Fuzzy recognition by the prokaryotic transcription factor HigA2 from Vibrio cholerae.

Disordered protein sequences can exhibit different binding modes, ranging from well-ordered folding-upon-binding to highly dynamic fuzzy binding. The primary function of the intrinsically disordered region of the antitoxin HigA2 from Vibrio cholerae is to neutralize HigB2 toxin through ultra-high-affinity folding-upon-binding interaction. Here, we show that the same intrinsically disordered region can also mediate fuzzy interactions with its operator DNA and, through interplay with the folded helix-turn-helix domain, regulates transcription from the higBA2 operon. NMR, SAXS, ITC and in vivo experiments converge towards a consistent picture where a specific set of residues in the intrinsically disordered region mediate electrostatic and hydrophobic interactions while "hovering" over the DNA operator. Sensitivity of the intrinsically disordered region to scrambling the sequence, position-specific contacts and absence of redundant, multivalent interactions, point towards a more specific type of fuzzy binding. Our work demonstrates how a bacterial regulator achieves dual functionality by utilizing two distinct interaction modes within the same disordered sequence.

Estimation of Peptide Helicity from Circular Dichroism Using the Ensemble Model.

An established method for the quantitation of the helix content in peptides using circular dichroism (CD) relies on the linear spectroscopic model. This model assumes an average value of the helix-length correction for all peptide conformers, irrespective of the length of the helical segment. Here we assess the validity of this approximation and introduce a more physically realistic ensemble-based analysis of the CD signal in which the length correction is assigned specifically to each ensemble conformer. We demonstrate that the linear model underestimates peptide helicity, with the difference depending on the ensemble composition. We developed a computer program that implements the ensemble model to estimate the peptide helicity. Using this model and the CD data set covering a broad range of helicities, we recalibrate CD baseline parameters and redetermine helix-coil parameters for the alanine-rich peptide. We show that the ensemble model leverages small differences in signal between conformers to extract more information from the experimental data, enabling the determination of several poorly defined quantities, such as the nucleation constant and heat capacity change associated with helix folding. Overall, the presented ensemble-based treatment of the CD signal, together with the recalibrated values of the spectroscopic baseline parameters, provides a coherent framework for the analysis of the peptide helix content.

Leucine Motifs Stabilize Residual Helical Structure in Disordered Proteins.

Many examples are known of regions of intrinsically disordered proteins that fold into α-helices upon binding to their targets. These helical binding motifs (HBMs) can be partially helical also in the unbound state, and this so-called residual structure can affect binding affinity and kinetics. To investigate the underlying mechanisms governing the formation of residual helical structure, we assembled a dataset of experimental helix contents of 65 peptides containing HBM that fold-upon-binding. The average residual helicity is 17% and increases to 60% upon target binding. The helix contents of residual and target-bound structures do not correlate, however the relative location of helix elements in both states shows a strong overlap. Compared to the general disordered regions, HBMs are enriched in amino acids with high helix preference and these residues are typically involved in target binding, explaining the overlap in helix positions. In particular, we find that leucine residues and leucine motifs in HBMs are the major contributors to helix stabilization and target-binding. For the two model peptides, we show that substitution of leucine motifs to other hydrophobic residues (valine or isoleucine) leads to reduction of residual helicity, supporting the role of leucine as helix stabilizer. From the three hydrophobic residues only leucine can efficiently stabilize residual helical structure. We suggest that the high occurrence of leucine motifs and a general preference for leucine at binding interfaces in HBMs can be explained by its unique ability to stabilize helical elements.

Binding-Induced Diversity of a Human Telomeric G-Quadruplex Stability Phase Space.

The structural polymorphism of G-quadruplex nucleic acids is an important factor in their recognition by proteins and small-molecule ligands. However, it is not clear why the binding of several ligands alters G-quadruplex topology. We addressed this question by following the (un)folding and binding of the human telomeric fragment 5'-(GGGTTA)3GGGT-3' (22GT) by calorimetry (DSC, ITC) and spectroscopy (CD). A thermodynamic analysis of the obtained data led to a detailed description of the topological phase space of stability (phase diagram) of 22GT and shows how it changes in the presence of a specific bisquinolinium ligand (360A). Various 1:1 and 2:1 ligand-quadruplex complexes were observed. With increasing temperature, the 1:1 complexes transformed into 2:1 complexes, which is attributed to the preferential binding of the ligand to the folding intermediates. Overall, the dissection of the thermodynamic parameters in combination with molecular modelling clarified the driving forces of the topological quadruplex transformations in a wide range of ligand concentrations and temperatures.

Analysis of Protein-DNA Interactions Using Isothermal Titration Calorimetry: Successes and Failures.

Isothermal titration calorimetry (ITC) is a golden standard for the characterization of protein-DNA binding affinities and allows direct assessment of the accompanying thermodynamic driving forces. Their interpretation can give insight into role of electrostatics, specificity of the DNA recognition, contribution of protein folding upon DNA binding and help to distinguish between minor and major groove binders. The main advantages of ITC are that the binding is measured in solution, and it requires no labeling of the samples, however, the method is not well suited for high-performance studies. Here we describe the sample preparation, a procedure to perform a typical ITC experiment, data analysis, and lastly discuss how to interpret the obtained thermodynamic parameters. In conclusion, we show examples of several unsuccessful ITC experiments and identify the underlying reasons for failed experiments. In most cases with a proper adjustment of the experimental setup, it was possible to obtain data appropriate for further analysis.

The free energy folding penalty accompanying binding of intrinsically disordered α-helical motifs.

Intrinsically disordered proteins (IDPs) are abundant in eukaryotic proteomes and preform critical roles in many cellular processes, most often through the association with globular proteins. Despite lacking a stable three-dimensional structure by themselves, they may acquire a defined conformation upon binding globular targets. The most common type of secondary structure acquired by these binding motifs entails formation of an α-helix. It has been hypothesized that such disorder-to-order transitions are associated with a significant free energy penalty due to IDP folding, which reduces the overall IDP-target affinity. However, the exact magnitude of IDP folding penalty in α-helical binding motifs has not been systematically estimated. Here, we report the folding penalty contributions for 30 IDPs undergoing folding-upon-binding and find that the average IDP folding penalty is +2.0 kcal/mol and ranges from 0.7 to 3.5 kcal/mol. We observe that the folding penalty scales approximately linearly with the change in IDP helicity upon binding, which provides a simple empirical way to estimate folding penalty. We analyze to what extent do pre-structuring and target-bound IDP dynamics (fuzziness) reduce the folding penalty and find that these effects combined, on average, reduce the folding cost by around half. Taken together, the presented analysis provides a quantitative basis for understanding the role of folding penalty in IDP-target interactions and introduces a method estimate this quantity. Estimation and reduction of IDP folding penalty may prove useful in the rational design of helix-stabilized inhibitors of IDP-target interactions. STATEMENT: The α-helical binding motifs are ubiquitous among the intrinsically disordered proteins (IDPs). Upon binding their targets, they undergo a disorder-to-order transition, which is accompanied by a significant folding penalty whose magnitude is generally not known. Here, we use recently developed statistical-thermodynamic model to estimate the folding penalties for 30 IDPs and clarify the roles of IDP pre-folding and bound-state dynamics in reducing the folding penalty.

Structural polymorphism of coiled-coils from the stalk domain of SARS-CoV-2 spike protein.

Spike trimer plays a key role in SARS-CoV-2 infection and vaccine development. It consists of a globular head and a flexible stalk domain that anchors the protein into the viral membrane. While the head domain has been extensively studied, the properties of the adjoining stalk are poorly understood. Here, we characterize the coiled-coil formation and thermodynamic stability of the stalk domain and its segments. We find that the N-terminal segment of the stalk does not form coiled-coils and remains disordered in solution. The C-terminal stalk segment forms a trimeric coiled-coil in solution, which becomes significantly stabilized in the context of the full-length stalk. Its crystal structure reveals a novel antiparallel tetramer coiled-coil with an unusual combination of a-d and e-a-d hydrophobic core packing. Structural analysis shows that a subset of hydrophobic residues stabilizes different coiled-coil structures: trimer, tetramer, and heterohexamer, underscoring a highly polymorphic nature of the SARS-CoV-2 stalk sequence.

Aggregation Time Machine: A Platform for the Prediction and Optimization of Long-Term Antibody Stability Using Short-Term Kinetic Analysis.

Monoclonal antibodies are the fastest growing class of therapeutics. However, aggregation limits their shelf life and can lead to adverse immune responses. Assessment and optimization of the long-term antibody stability are therefore key challenges in the biologic drug development. Here, we present a platform based on the analysis of temperature-dependent aggregation data that can dramatically shorten the assessment of the long-term aggregation stability and thus accelerate the optimization of antibody formulations. For a set of antibodies used in the therapeutic areas from oncology to rheumatology and osteoporosis, we obtain an accurate prediction of aggregate fractions for up to three years using the data obtained on a much shorter time scale. Significantly, the strategy combining kinetic and thermodynamic analysis not only contributes to a better understanding of the molecular mechanisms of antibody aggregation but has already proven to be very effective in the development and production of biological therapeutics.

The sequence-ensemble relationship in fuzzy protein complexes.

Intrinsically disordered proteins (IDPs) interact with globular proteins through a variety of mechanisms, resulting in the structurally heterogeneous ensembles known as fuzzy complexes. While there exists a reasonable comprehension on how IDP sequence determines the unbound IDP ensemble, little is known about what shapes the structural characteristics of IDPs bound to their targets. Using a statistical thermodynamic model, we show that the target-bound ensembles are determined by a simple code that combines the IDP sequence and the distribution of IDP-target interaction hotspots. These two parameters define the conformational space of target-bound IDPs and rationalize the observed structural heterogeneity of fuzzy complexes. The presented model successfully reproduces the dynamical signatures of target-bound IDPs from the NMR relaxation experiments as well as the changes of interaction affinity and the IDP helicity induced by mutations. The model explains how the target-bound IDP ensemble adapts to mutations in order to achieve an optimal balance between conformational freedom and interaction energy. Taken together, the presented sequence-ensemble relationship of fuzzy complexes explains the different manifestations of IDP disorder in folding-upon-binding processes.

Unraveling the Thermodynamics of Ultra-tight Binding of Intrinsically Disordered Proteins.

Protein interactions mediated by the intrinsically disordered proteins (IDPs) are generally associated with lower affinities compared to those between globular proteins. Here, we characterize the association between the intrinsically disordered HigA2 antitoxin and its globular target HigB2 toxin from Vibrio cholerae using competition ITC experiments. We demonstrate that this interaction reaches one of the highest affinities reported for IDP-target systems (K D = 3 pM) and can be entirely attributed to a short, 20-residue-long interaction motif that folds into α-helix upon binding. We perform an experimentally based decomposition of the IDP-target association parameters into folding and binding contributions, which allows a direct comparison of the binding contribution with those from globular ultra-high affinity binders. We find that the HigA2-HigB2 interface is energy optimized to a similar extent as the interfaces of globular ultra-high affinity complexes, such as barnase-barstar. Evaluation of other ultra-high affinity IDP-target systems shows that a strategy based on entropy optimization can also achieve comparably high, picomolar affinities. Taken together, these examples show how IDP-target interactions achieve picomolar affinities either through enthalpy optimization (HigA2-HigB2), resembling the ultra-high affinity binding of globular proteins, or via bound-state fuzziness and entropy optimization (CcdA-CcdB, histone H1-prothymosin α).

Conditional Cooperativity in DNA Minor-Groove Recognition by Oligopeptides.

The recognition of specific DNA sequences in processes such as transcription is associated with a cooperative binding of proteins. Some transcription regulation mechanisms involve additional proteins that can influence the binding cooperativity by acting as corepressors or coactivators. In a conditional cooperativity mechanism, the same protein can induce binding cooperativity at one concentration and inhibit it at another. Here, we use calorimetric (ITC) and spectroscopic (UV, CD) experiments to show that such conditional cooperativity can also be achieved by the small DNA-directed oligopeptides distamycin and netropsin. Using a global thermodynamic analysis of the observed binding and (un)folding processes, we calculate the phase diagrams for this system, which show that distamycin binding cooperativity is more pronounced at lower temperatures and can be first induced and then reduced by increasing the netropsin or/and Na+ ion concentration. A molecular interpretation of this phenomenon is suggested.

A nanobody toolbox targeting dimeric coiled-coil modules for functionalization of designed protein origami structures.

Coiled-coil (CC) dimers are widely used in protein design because of their modularity and well-understood sequence-structure relationship. In CC protein origami design, a polypeptide chain is assembled from a defined sequence of CC building segments that determine the self-assembly of protein cages into polyhedral shapes, such as the tetrahedron, triangular prism, or four-sided pyramid. However, a targeted functionalization of the CC modules could significantly expand the versatility of protein origami scaffolds. Here, we describe a panel of single-chain camelid antibodies (nanobodies) directed against different CC modules of a de novo designed protein origami tetrahedron. We show that these nanobodies are able to recognize the same CC modules in different polyhedral contexts, such as isolated CC dimers, tetrahedra, triangular prisms, or trigonal bipyramids, thereby extending the ability to functionalize polyhedra with nanobodies in a desired stoichiometry. Crystal structures of five nanobody-CC complexes in combination with small-angle X-ray scattering show binding interactions between nanobodies and CC dimers forming the edges of a tetrahedron with the nanobody entering the tetrahedral cavity. Furthermore, we identified a pair of allosteric nanobodies in which the binding to the distant epitopes on the antiparallel homodimeric APH CC is coupled via a strong positive cooperativity. A toolbox of well-characterized nanobodies specific for CC modules provides a unique tool to target defined sites in the designed protein structures, thus opening numerous opportunities for the functionalization of CC protein origami polyhedra or CC-based bionanomaterials.

G-Quadruplex Stability from DSC Measurements.

Guanine-rich DNA oligonucleotides can adopt G-quadruplex (G4) structures in the presence of specific cations. Folding and unfolding of G4 can be characterized thermodynamically, providing the information on the stability of various G4 conformations. We show how the relevant thermodynamic and sometimes kinetic parameters are obtained by employing differential scanning calorimetry (DSC) and global fitting of an appropriate model to the DSC data. We demonstrate that best-fit values of the thermodynamic parameters can be interpreted in terms of the driving forces accompanying the G4 folding/interconversion and how they are translated into the phase diagrams, which provide an elegant description of the G4 phase space over a wide range of solution conditions.

Thermodynamic Stability of the Transcription Regulator PaaR2 from Escherichia coli O157:H7.

PaaR2 is a putative transcription regulator encoded by a three-component parDE-like toxin-antitoxin module from Escherichia coli O157:H7. Although this module's toxin, antitoxin, and toxin-antitoxin complex have been more thoroughly investigated, little remains known about its transcription regulator PaaR2. Using a wide range of biophysical techniques (circular dichroism spectroscopy, size-exclusion chromatography-multiangle laser light scattering, dynamic light scattering, small-angle x-ray scattering, and native mass spectrometry), we demonstrate that PaaR2 mainly consists of α-helices and displays a concentration-dependent octameric build-up in solution and that this octamer contains a global shape that is significantly nonspherical. Thermal unfolding of PaaR2 is reversible and displays several transitions, suggesting a complex unfolding mechanism. The unfolding data obtained from spectroscopic and calorimetric methods were combined into a unifying thermodynamic model, which suggests a five-state unfolding trajectory. Furthermore, the model allows the calculation of a stability phase diagram, which shows that, under physiological conditions, PaaR2 mainly exists as a dimer that can swiftly oligomerize into an octamer depending on local protein concentrations. These findings, based on a thorough biophysical and thermodynamic analysis of PaaR2, may provide important insights into biological function such as DNA binding and transcriptional regulation.

Energetic Basis of AGCGA-Rich DNA Folding into a Tetrahelical Structure.

It was recently discovered that, besides well-known G-quadruplexes and i-motifs, DNA may adopt another type of noncanonical structure called AGCGA-quadruplexes. Here, the folding of the VK2 fragment from the regulatory region of the PLEKHG3 gene is studied and, for the first time, the energetic contributions that stabilize this unique fold are described. Similarly to the B-DNA, it is stabilized by hydrophobic desolvation and, in contrast to G-quadruplexes, also by specific binding of water molecules. Compared to B-DNA, VK2 folding is enthalpically less favorable due to poorer base-stacking interactions, resulting in substantial conformational flexibility. This entropically favorable conformational "breathing" stabilizes the AGCGA-quadruplexes. In conclusion, AGCGA-quadruplexes have a distinguishing thermodynamic fingerprint and the corresponding driving forces enabling their folding are consistent with the observed structural features.

Structural Basis of Epitope Recognition by Heavy-Chain Camelid Antibodies.

Truncated versions of heavy-chain antibodies (HCAbs) from camelids, also termed nanobodies, comprise only one-tenth the mass of conventional antibodies, yet retain similar, high binding affinities for the antigens. Here we analyze a large data set of nanobody-antigen crystal structures and investigate how nanobody-antigen recognition compares to the one by conventional antibodies. We find that nanobody paratopes are enriched in aromatic residues just like conventional antibodies, but additionally, they also bear a more hydrophobic character. Most striking differences were observed in the characteristics of the antigen's epitope. Unlike conventional antibodies, nanobodies bind to more rigid, concave, conserved and structured epitopes enriched with aromatic residues. Nanobodies establish fewer interactions with the antigens compared to conventional antibodies, and we speculate that high binding affinities are achieved due to less unfavorable conformational and more favorable solvation entropy contributions. We observed that interactions with antigen are mediated not only by three CDR loops but also by numerous residues from the nanobody framework. These residues are not distributed uniformly; rather, they are concentrated into four structurally distinct regions and mediate mostly charged interactions. Our findings suggest that in some respects nanobody-antigen interactions are more similar to the general protein-protein interactions rather than antibody-antigen interactions.

The Thermodynamic Basis of the Fuzzy Interaction of an Intrinsically Disordered Protein.

Many intrinsically disordered proteins (IDP) that fold upon binding retain conformational heterogeneity in IDP-target complexes. The thermodynamics of such fuzzy interactions is poorly understood. Herein we introduce a thermodynamic framework, based on analysis of ITC and CD spectroscopy data, that provides experimental descriptions of IDP association in terms of folding and binding contributions which can be predicted using sequence folding propensities and molecular modeling. We show how IDP can modulate the entropy and enthalpy by adapting their bound-state structural ensemble to achieve optimal binding. This is explained in terms of a free-energy landscape that provides the relationship between free-energy, sequence folding propensity, and disorder. The observed "fuzzy" behavior is possible because of IDP flexibility and also because backbone and side-chain interactions are, to some extent, energetically decoupled allowing IDP to minimize energetically unfavorable folding.

Molecular mechanism governing ratio-dependent transcription regulation in the ccdAB operon.

Bacteria can become transiently tolerant to several classes of antibiotics. This phenomenon known as persistence is regulated by small genetic elements called toxin-antitoxin modules with intricate yet often poorly understood self-regulatory features. Here, we describe the structures of molecular complexes and interactions that drive the transcription regulation of the ccdAB toxin-antitoxin module. Low specificity and affinity of the antitoxin CcdA2 for individual binding sites on the operator are enhanced by the toxin CcdB2, which bridges the CcdA2 dimers. This results in a unique extended repressing complex that spirals around the operator and presents equally spaced DNA binding sites. The multivalency of binding sites induces a digital on-off switch for transcription, regulated by the toxin:antitoxin ratio. The ratio at which this switch occurs is modulated by non-specific interactions with the excess chromosomal DNA. Altogether, we present the molecular mechanisms underlying the ratio-dependent transcriptional regulation of the ccdAB operon.

Ribosome-dependent Vibrio cholerae mRNAse HigB2 is regulated by a ß-strand sliding mechanism.

Toxin-antitoxin (TA) modules are small operons involved in bacterial stress response and persistence. higBA operons form a family of TA modules with an inverted gene organization and a toxin belonging to the RelE/ParE superfamily. Here, we present the crystal structures of chromosomally encoded Vibrio cholerae antitoxin (VcHigA2), toxin (VcHigB2) and their complex, which show significant differences in structure and mechanisms of function compared to the higBA module from plasmid Rts1, the defining member of the family. The VcHigB2 is more closely related to Escherichia coli RelE both in terms of overall structure and the organization of its active site. VcHigB2 is neutralized by VcHigA2, a modular protein with an N-terminal intrinsically disordered toxin-neutralizing segment followed by a C-terminal helix-turn-helix dimerization and DNA binding domain. VcHigA2 binds VcHigB2 with picomolar affinity, which is mainly a consequence of entropically favorable de-solvation of a large hydrophobic binding interface and enthalpically favorable folding of the N-terminal domain into an α-helix followed by a ß-strand. This interaction displaces helix α3 of VcHigB2 and at the same time induces a one-residue shift in the register of ß-strand ß3, thereby flipping the catalytically important Arg64 out of the active site.

Unraveling the Thermodynamics of the Folding and Interconversion of Human Telomere G-Quadruplexes.

Why human telomere DNA fragments fold into different G-quadruplex structures with parallel, hybrid, and antiparallel strand orientations depending on the temperature and concentration of co-solutes remains poorly understood. Similarly, the formation of intermediate structures along the folding or interconversion pathways is not well understood. Herein, we address these questions by introducing a conceptual framework, based on the global thermodynamic analysis of DSC and CD spectroscopy data, which led to a detailed description of the topological phase space (phase diagram) of the stability of the human telomere fragment 5'-AGGG(TTAGGG)3 -3' (Tel22). This framework clarifies the driving forces of quadruplex folding and interconversion processes over a wide range of temperatures and ion (K(+) , Na(+) ) and polyethylene glycol (PEG) concentrations and demonstrates their linkage to the human telomere DNA structural features.